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' Abstract 

We have developed a Monte Carlo event generator for non-resonant diphoton (77) production 
^-j, at hadron collisions in the framework of GR@PPA, which consistently includes processes having 

additional one jet radiation. The possible double count problem in the generation of radiative 
processes is avoided by using the LLL subtraction method that we have applied to the weak- 
boson production processes. The subtraction method has been extended to the final-state QED 
divergence that appears in the qg — > 77 + q process. Because a parton shower (PS) which 
regularizes the subtracted QED divergence is still under development, we tried to use PYTHIA for 
the generation of the fragmentation events to restore the subtracted components. The simulation 
employing the "old" PS of PYTHIA shows a reasonable matching with the GR@PPA events, and 
O |. the combined event sample shows a result in reasonable agreement with ResBos. We found that 

the contribution from qg — > 77 + q is significant in the LHC condition. This event generator must 
be useful for the background studies in low-mass Higgs boson searches at LHC. 
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1 Introduction 

> " 

Diphoton (77) production is one of most promising channels for the discovery of the Higgs boson 
having a relatively small invariant mass (< 130 GeV/c 2 ) at LHC. However, the measurement must 
ly-j j suffer from a large irreducible diphoton background produced via non-resonant QED interactions. 

It is necessary to understand the properties of this background, not only for the discovery of the 
Higgs boson but also for detailed studies after the discovery. The identification of photons is rather 
complicated and largely dependent on the detector performance. For instance, a certain isolation 
condition has to be required in order to reduce large contamination of ir° from hadron jets. The 
performance of this selection depends on the details of detector responses and is hard to be evaluated 
analytically. Therefore, it is strongly desired to provide theoretical predictions in the form of Monte 
^ . Carlo (MC) event generators. 

The lowest order process for non-resonant diphoton production is very simple as shown in Fig. [T] 
(a). Despite that, the next-to-leading order (NLO) correction to this process is known to be very large 
[1] [2] . The large correction is predominantly due to the contribution from real radiation processes 
illustrated in Fig. Q] (b) to (d). While the contribution from gluon radiation processes shown in Fig. 
[1] (b) is not very significant, those from quark radiation processes shown in Fig. [1] (c) and (d) may 
become large due to a very large gluon density inside protons. It must be necessary to include these 
processes in order to provide a realistic simulation. 

In this report, we describe the MC event generator for non-resonant diphoton production at hadron 
collisions that we have developed. The program has been developed in the framework of the GR@PPA 
event generator [3], and supports the generation of radiative processes in Fig.[T](b) to (d) together with 
the lowest-order process in Fig.[T](a). Though the generator includes radiative processes, it is not fully 
including NLO corrections. Non-divergent terms in soft/virtual corrections are yet to be included. In 
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Figure 1: Typical Feynman diagrams for non- resonant diphoton production at hadron collisions: (a) 
the lowest order, (b) a gluon radiation process, (c) a quark radiation process, and (d) another quark 
radiation process. The processes (b) and (c) have initial-state QCD divergences, while (d) has a 
final-state QED divergence. 

any case, the radiative processes have various divergences which we need to regularize. The initial- 
state QCD divergences can be regularized using the method that we have applied to weak-boson 
production processes, where divergent terms are numerically subtracted from the matrix elements 
of radiative processes (the LLL subtraction) [4j \5\ [6] . The subtracted contributions are restored by 
combining with non-radiative processes to which a parton shower (PS) is applied. We can avoid the 
double count problem by the subtraction and naturally regularize the divergences as a result of the 
multiple radiation in PS. 

The quark radiation process illustrated in Fig. Q] (c) and (d) not only has an initial-state QCD 
divergence but also has a QED divergence in the final state. We have extended the method applied 
to the initial-state QCD to this final-state QED divergence. The extension of the subtraction is 
straightforward, while the preparation of an appropriate PS is not trivial since it has to support 
QED radiations together with QCD. This PS has to be applied to the quark in the final state of the 
9.9 ~> 7 + 9 process to radiate photons from the final-state quark based on a collinear approximation 
(the fragmentation process). It is desired to have a capability to force a hard photon radiation in 
this PS since we are interested in those events having two hard photons in the final state. Such a 
final-state PS is still in development. Instead, we try to use the PYTHIA PS [7] for simulating the 
fragmentation process in this report. Though the PYTHIA PS is capable of radiating photons, it 
does not have a mechanism to force a hard photon radiation. We need to repeat the generation until 
we observe a sufficiently hard photon in the generated event. We use the so-called "old" PS in the 
present study. The "new" PS is not used because we still have some questions on its behavior. 

We sometimes find reports in which the detection efficiency and acceptance for the diphoton 
measurements are evaluated using event generators for the lowest-order process and fragmentation 
processes. Such evaluations are not self-consistent. Parton showers used in the simulation of fragmen- 
tation processes have a certain energy scale defining the maximum hardness of the parton radiation. 
The results depend on this arbitrarily chosen energy scale since those radiations exceeding this scale 
are ignored. Non-collinear contributions are also ignored. These ignored contributions may be small 
in many processes compared to the contributions taken into consideration. However, as we will show 
in this report, they become comparable to the lowest-order contribution in the diphoton production. 
It is necessary to include a simulation based on the exact matrix elements for the qg — ¥ 77 + 9 process 
in order to make a reliable evaluation. 

The fragmentation process that we take into consideration in this report is the so-called "single" 
fragmentation. The "double" fragmentation in which two photons are radiated from final-state par- 
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tons, for instance, from quarks in the gg — > qq process are not supported. We need to introduce "77 + 
2 jets" production processes in order to construct a consistent event generator including the "double" 
fragmentation. Besides, gg — ¥ 77 and its higher orders are not included at present. 

We require a typical kinematical condition for the Higgs-boson search at LHC through the present 
study because we are interested in its background; that is 

Pt{h) > 40 GeV/c, p T (7a) > 25 GeV/c, 
|77(7)| < 2.5, AR{jj) > 0.4 

80 < m 77 < 140 GeV/c 2 . (1) 

We apply an asymmetric requirement to the transverse momenta (pr) of photons with respect to 
the incident beam direction. The requirement on the pseudorapidity (77) is common to the two 
photons. In addition, though this is not effective for real diphoton events now we consider, we require 
a sufficient AR separation between the two photons, where AR is defined from the differences in the 
pseudorapidity (7/) and the azimuthal angle (</>) as AR 2 = Aij 2 + A(f> 2 . Finally, the invariant mass 
of the two photons (77t 77 ) is restricted within the range that we are interested in. These conditions 
are required after completing the simulation down to the hadron level. Looser constraints are applied 
at the event generation in order to avoid any bias. The simulations are carried out for the design 
condition of LHC, proton-proton collisions at a center-of-mass (cm) energy of 14 TeV. 

This report is organized as follows: the extension of the limited leading-log (LLL) subtraction 
method to the final-state QED divergence is described in Section 2, and the simulation of the frag- 
mentation process employing the PYTHIA PS is described in Section 3. A typical isolation cut is 
simulated in Section 4. The results from the combined event simulation is presented in Section 5, and 
the discussions are concluded in Section 6. 



2 Final-state QED LLL subtraction 



We approximate the final-state QED divergence in the qg — > 77 + q process as 
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where e q is equal to 1/3 for down- type quarks and 2/3 for up- type quarks. We evaluate Eq. ^ for 
two possible combinations of the quark and a photon in the final state, and numerically subtract 
them from the exact matrix element of qg — > 77 + q, together with the LLL term for the initial-state 
QCD divergence. The parameters Q 2 and z are defined by using the sum of the energy (E q7 ) and the 
momentum (p 97 ) of the considered g-7 pair. They are defined in the cm frame of the qg — > 77 + q 
event as 



Q = E 2 



Pa 



and 



PL {E ql + Pql )+Q 2 /2 



(5) 



(6) 



P qi {E qi +p qi )+Q 2 

where pl is the momentum component of the q in parallel to the summed momentum. Equation 
© is defined so that z should represent the momentum fraction at the infinite-momentum limit, 
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Equation ([3]) is slightly different from the radiation function for the initial state [4] ; z has vanished 
in the denominator. This is because we assume that the cm energy does not change before and after 
the radiation. The mapping to the non-radiative process, qg — >• 7 + q, which is necessary to perform 
in order to evaluate the LLL term in Eq. ^ , is also defined according to this assumption. We define 
the momentum of the quark in the final state of the mapped qg — > 7 + q event by the sum of the 
momenta of the g-7 pair in consideration. Since all particles in the qg — > 7 + q event are assumed to be 
on-shell, the invariant mass of the final state evaluated from thus defined momentum becomes smaller 
than that of the initial state. In order to compensate for this decrease, we increase the overall scale 
of the momenta of the final-state particles. The Q 2 value defined in Eq. ([5]) is also increased with the 
same scaling factor. The z value is independent of this rescaling. These details, the definitions in Eqs. 
([5]) and © and the subsequent momentum adjustment, are not universal. They have been chosen in 
order to achieve a good matching with the parton shower (PS) that we are developing. Therefore, the 
application of the PYTHIA PS may result in a certain mismatch, though the effect of such details 
must vanish at the limit where the radiative cross section diverges. 

An event generator implementing the above subtraction has been developed in the framework of 
GR@PPA, and the generation was tried for proton-proton collisions with a cm energy of 14 TeV with 
CTEQ6L1 [8] used for PDF. The energy scale was defined as 

V 2 = \pt(h) ~pt(72)| 2 /4, (7) 

where prdi) denotes the transverse momentum vector of the two photons. This definition is equivalent 
to the ordinary definition, fi = pt{i), for qq — > 77. We used the same definition for the renormalization 
and factorization scales. The energy scale for the initial-state PS must be equal to the factorization 
scale in our method. Though it is not necessary, we adopted the same definition for the final-state 
PS. Thus, all the energy scales used in the event generation were identical in this test. 

The LLL subtraction is limited by the 9 function in Eq. @. This is because the implementation 
of the PS to be applied to the non-radiative process is limited by a certain energy scale. Thus, the 
energy scale used in Eq. ([2]), {J-fsr, must be equal to the one defined for the mapped qg — > 7 + q event 
in order to achieve a good matching. In the present study, we define hfsr to be equal to the px of 
the mapped qg — > 7 + q event. 

One of the simulation results is shown in Fig. [3] In this figure, we have plotted the distribution 
of AR between the photon and the quark in the final state of the qg — > 77 + q events, where q 
represents any quark or anti-quark up to the b(b) quark. We obtain two values since there are two 
photons in the final state. We take the smaller one as the AR(j-q). Parton showers are yet to be 
applied but the initial-state QCD LLL subtraction is already applied in this simulation. The selection 
condition in Eq. ([1]) is applied to the generated events. The distribution directly derived from the 
matrix element before applying the final-state LLL subtraction is shown with a solid histogram in 
the figure for comparison. We applied a cut of AR(j-q) > 0.2 since this distribution is divergent at 
AR(j-q) = 0. 

Since the LLL subtraction is unphysical, we obtain negative-weight events as well as ordinary 
positive- weight events when we apply the subtraction. The event weights are always equal to +1 or 
— 1 because BASES/SPRING [HI [TU] automatically unweights the events. Therefore, we can obtain 
the desired distribution by subtracting the number of negative-weight events from that of positive- 
weight events in each histogram bin. The open circles (triangles) in the figure show the distribution of 
positive(negative)-weight events, and the final distribution is shown with filled circles. We can see that 
the subtraction is effective only in a relatively small AR("/-q) region as expected. The distribution 
at large AR(^-q) is not altered by the subtraction. We can also see that the distributions after the 
subtraction converge to finite values as AR(j-q) — »■ 0; not only the final result converges to zero, but 
also the positive and negative weight events converge to a finite value at this limit. These facts show 
that the subtraction is done properly at least near the divergent limit. We have applied a small cutoff, 
AR(j-q) > 0.01, in the event generation for numerical stability. We can see that the effect of this 
cutoff is negligible. 
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Figure 2: AR^j-q) distribution of the simulated qg — > 77 + q events satisfying the kinematical con- 
dition in Eq. (JlJ. The solid histogram shows the distribution before applying the final-state QED 
LLL subtraction, where AR(j-q) > 0.2 is required to cutoff the divergence. Open circles show the 
distribution of positive- weight events after the subtraction, and open triangles show that of negative- 
weight events. The final result after applying the final-state QED LLL subtraction is shown with filled 
circles. The initial-state QCD LLL subtraction is already applied while parton showers are yet to be 
applied in this simulation. 

3 Fragmentation process by PYTHIA PS 

The subtracted LLL terms must be restored by non-radiative processes to which an appropriate parton 
shower (PS) is applied. Because our PS is still under development, we try to use the PYTHIA PS [7] 
in the present study. We used the so-called "old" PS because this model looks similar to the PS that 
we are developing. We used PYTHIA 6.423 for this study. We generated qq — >• 37 and qg — > qj events 
by setting as msel = 0,msub(14) = 1, and msub(29) = 1 in the LHC condition at the design energy, 
14 TeV. The qq — > 37 interaction was turned on for completeness. Though it is very inefficient to find 
hard photons in the gluon fragmentation, turning on this process is harmless since the production 
cross section is small compared to qg — > qj. 

For PDF, CTEQ6L1 was applied with the help of LHAPDF/LHAGLUE [TT] in the LHAPDF 
5.8.4 distribution by setting as mstp(51) = 10042 and mstp(52) = 2 in PYTHIA. The final-state 
particles were allowed to produce in the region \r)\ < 3.0 with the transverse momenta of pr > 30 
GeV/c at the on-shell parton level by appropriately setting the ckin parameters. Though this px cut 
may look tight compared to the condition in Eq. ([T]), it is safe enough because the lower-pT photons 
are predominantly produced in the fragmentation. Since we do not want to change the preset energy 
scales, we set the scale parameters as parp(67) = 1.0 and parp(71) = 1.0. All the energy scales 
must have been set equal to the pt of the on-shell parton level interaction with this setting. In 
addition, for safety, we explicitly disabled the matrix-element correction by setting as mstp(68) = 0. 
The other parameters were left unchanged so that the default "old" PS and further simulations down 
to the hadron level should be applied according to the default setting. 

The event generation with the above setting was repeated and hard photons were looked for in the 
fragmentation of the final-state partons. The probability to find hard photons in PS is very small. 
The efficiency to find events satisfying the condition in Eq. ([T]) was 2.5 x 10~ 4 . In order to improve 
the efficiency, we enhanced the QED radiation by a factor of 10 by setting as mst j (41) = 10 and 
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Figure 3: Ai2(7-jet) distribution of the simulated 77+ jet events satisfying the kinematical condition 
in Eq. (jTJ) , where "jet" represents the quark or gluon in the final state. The solid histogram shows 
the distribution of the fragmentation events simulated by PYTHIA. The dashed histogram shows the 
distribution of the qg — » 77 + q events generated by GR@PPA in which the initial-state and final- 
state LLL subtractions are fully applied. In both simulations, the event selection has been applied 
to those events fully simulated down to the hadron level, while Ai?(7-jet) has been derived from the 
information of the reconstructed or original 77+ jet events at the on-shell parton level. The sum of the 
two distributions is plotted with hllcd circles. The dotted histogram shows the distribution directly 
derived from the qg — > 77 + q matrix element before applying the final-state QED LLL subtraction 
(same as the solid histogram in Fig. [5]). The distribution of the qq — > 77 + g events is shown with a 
dot-dashed histogram for comparison. 

parj (84) = 10.0. The obtained results were corrected for this enhancement. It is not recommended 
in the manual to apply larger enhancement factors because the multiple photon radiation effect may 
become significant. By the way, it should be noted that the efficiency is still at the level of 10 -3 even 
with this enhancement. 

In order to investigate the matching with the subtracted GR@PPA simulation results, it is better 
to reconstruct the corresponding qg — > 77 + q events at the on-shell parton level from full PYTHIA 
simulation events in which a hard fragmentation photon is observed. This reconstruction can be done 
unambiguously by using the information of the initial-state partons initiated the hard interaction, 
without using the information of the final-state partons or hadrons. The reconstruction using the 
final-state partons or hadrons is dangerous because their origin is ambiguous. The initial-state partons 
initiated the hard interaction can be found in the documentation lines of the PYTHIA event record 
as partons having no child. We first boost the prompt and fragmentation photons to the cm frame 
of these partons in which the initial-state parton momenta are aligned along the beam direction. 
We determine the momentum of the remnant on-shell parton from the photon momenta to balance 
the total momentum. The flavor of the remnant parton can be determined from the flavor of the 
initial-state partons. The total energy of the final state derived from thus determined momenta is 
usually different from the initial-state cm energy as a result of the application of PS and hadronization 
simulations. We adjust the overall scale of the momenta of the final-state particles to match the total 
energy as is done in the mapping to non-radiative events in the LLL subtraction. We then boost the 
reconstructed parton and photon momenta to the laboratory frame by using the momentum fraction 
information in pari (33) and pari (34) originally used for the generation of the hard interaction. 
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The Ai?(7-jet) distribution of the simulated fragmentation events is shown with a solid histogram 
in Fig. [3J in which the selection condition in Eq. (JT|) was applied to the simulated hadron-level events, 
while Ai?(7-jet) was evaluated using the reconstructed on-shell parton level qg — > 77 + q information. 
In order to compare with this result, we applied a full simulation down to the hadron level to the 
GR@PPA simulation events described in the previous section. The events were generated with a 
looser kinematical constraint, pt{ii) > 30 GeV/c,pT(72) > 20 GeV/c, and 1 77 ('y) I — 3.0, and the 
initial-state and final-state parton showers were applied down to Qq — (4.6 GeV) 2 in GRQPPA. 
The generated events were passed to PYTHIA in order to simulate parton showers at smaller Q 2 
and hadronization/decays. The default setting was unchanged in the PYTHIA simulation, except 
for the setting of parp( 67) = 1 . and parp(71) = 1.0. The selection condition in Eq. (TTJ) was then 
applied to the simulated hadron-level events. The 7-jet separation, Ai?(7-jet), was evaluated from the 
original on-shell parton level information for the selected events. The obtained Ai?(7-jet) distribution 
is presented with a dashed histogram in Fig. [3] 

The obtained two distributions seem to be smoothly connected with some overlap around the 
boundary. We expect that the sum of them should reproduce the spectrum before the final-state 
subtraction shown with a dotted histogram around the boundary, since we have chosen the boundary, 
Q 2 — n 2 F: in a hard radiation region. The sum should gradually get apart from the prediction without 
the subtraction to approach the prediction from the fragmentation as Ai?(7-jet) becomes smaller. 
The sum shown with closed circles in Fig. [3] behaves nearly as expected. However, we can see a small 
dip at the boundary. This dip must be due to a certain mismatch in the hard radiation kinematics. 
Furthermore, the distribution of the fragmentation events seems to be a little bit smaller than we 
expect. The distribution looks better if we enhance the fragmentation by a factor of about 20%. We 
expect that we can achieve a better matching with the parton shower that we are developing. 

People may worry that the drop in the first bin of the fragmentation event distribution in Fig. [3] 
looks unnatural. This drop is caused by a cutoff of the QED PS in PYTHIA. The cutoff is set to 1 GeV 
in terms of Q by the default, to be equal to that of the QCD PS. Possible QCD phenomena at smaller 
Q values are simulated otherwise, for instance, by the multiple interaction, while no such simulations 
are implemented for QED. Thus, there is no reason to stop the QED PS at Q = 1 GeV. Hard 
photons can be radiated from branches having smaller Q values since the momentum determination is 
independent of the choice of Q. Actually, if we decrease this cutoff to 0.3 GeV by setting as par j (83) 
= . 3 in PYTHIA, the concentration in the first bin increases significantly. Though the concentration 
will further increase if we apply smaller cutoff values, there must be a certain limitation due to non- 
perturbative effects of QCD. The picture in which a quark branches to a quark and a photon should 
break at Q values smaller than the QCD cutoff. The treatment of such small-Q 2 branches is one of 
the subjects in our development. 

In Fig. [3j the contribution from qq — > 77 + g events is shown with a dot-dashed histogram for 
comparison. The initial-state QCD LLL subtraction is applied also to this process. This process does 
not have any final-state divergence. We can see that the contribution from this process is always 
smaller than qg — > 77 + q by nearly a factor of five. 

4 Isolation cut 

It is necessary in real experiments to require a certain isolation condition in the identification of 
photons, in order to reduce the huge background from ir° in hadron jets. It is difficult to reproduce 
the cuts applied by experiments with parton-level simulations. This is the main reason why hadron- 
level event generators are desired. If hadron-level events can be generated consistently, the generated 
events can be passed to detector simulations for more detailed studies including detector responses. 

In this section we try to simulate a typical isolation condition using the events simulated down to 
the hadron level. The condition that we require is defined by using a cone Et defined as 

-Et.coiic = ^ Et, (8) 

AR<R iBO wrt. 7 

where Et is the transverse component of the energy of particles with respect to the beam direction, 
and the sum is taken over all particles inside a given AR cone around the photon, excluding the 
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Figure 4: Ai?(7-jet) distribution of the simulated 77 + jet events. The isolation cut described in the 
text is applied to the events from which the distributions in Fig. [3] are derived. The notation of the 
histograms and the plot is the same as Fig. [3] 



photon under the study and neutrinos. We take as i?i SO = 0.4 and require that i?T,cono < 15 GeV. 
We applied this cut to the events from which the distributions in Fig. [3] were derived. The Ai?(7-jet) 
distribution after applying the cut is shown in Fig. [4j We can see that the isolation cut strongly 
suppresses the fragmentation events in Ai?(7-jet) < 0.4, as expected, while it has a very small impact 
to the LLL subtracted qg —> 77 + q and qq — > 77 + g events. 

Figure [3] shows the px distribution of the "jets" in the on-shell parton level 77 + jet events re- 
constructed for the fragmentation events in the region, Ai?(7-jet) < 0.4. The distribution after the 
cut should sharply drop at 15 GeV/c if the isolation cut at the on-shell parton level could perfectly 
reproduce the cut applied to the hadron-level events. We can see an apparent smearing around the 
ideal cut value. This smearing would result in an inaccuracy of the approximation in analytical eval- 
uations, such as DIPHOX and ResBos. In any case, the fact that the distribution drops around the 
ideal cut value, as we expect, must imply that the "reconstruction" of the on-shell parton level events 
is reasonably done. 



5 Combined simulation 



We can obtain a consistent simulation sample by combining the fragmentation events with the LLL 
subtracted qg — > 77 + q and 99—5-77 + 5 events, together with the lowest-order qq — > 77 events. 
The lowest-order events were simultaneously generated by GR@PPA with the LLL subtracted events, 
to restore the subtracted initial-state divergent components. A loose kinematic condition, pt(ii) > 
30 GeV/c,j?T(72) > 20 GeV/c, and (77(7)1 ^ 3-0, was commonly applied in the on-shell parton level 
event generation in GR@PPA. The event generation was done for the LHC design condition, pp 
collisions at 14 TeV, with the energy scale definition in Eq. ([7|). Built-in parton shower simulations 
were applied in GR@PPA down to Qq = (4.6 GeV) 2 . Further small-Q 2 phenomena were simulated 
by PYTHIA 6.423, with its default setting but parp (67) = 1 . and parp(71) = 1 . 0. The selection 
condition in Eq. (TTJ) was applied to thus simulated hadron-level events. The fragmentation events 
were separately simulated with PYTHIA as described in previous sections. Finally, the isolation cut 
described in the previous section was applied to the simulated events. 
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Figure 5: pr distribution the "jets" in the on-shell parton level 77 + jet events reconstructed for 
the fragmentation events in Ai?(7-jet) < 0.4. The solid histogram shows the distribution before the 
isolation cut, and the dashed histogram shows the distribution after the cut, i?T.conc < 15 GeV. 

Figure [5] shows the pr distribution of the diphoton (77) system in the simulated events. The 
obtained distribution is compared with the prediction from RcsBos 2 . ResBos provides us with an 
NLO prediction in which soft QCD radiations are resummed. It is considered to be most reliable 
at present as far as the pt(u) distribution is concerned. Unfortunately, since the resummation can 
predict inclusive properties only, ResBos cannot provide exclusive event information. It should be 
noted that the contribution from the gg — > 77 process is not included in the ResBos prediction since 
it is yet to be included in our simulation. The ResBos prediction shown in the figure is smaller than the 
result presented in the original paper because the to 77 constraint in Eq. ([lj is additionally required. 
We have confirmed that we can obtain a prediction close to the original result if the m 77 constraint 
is excluded. 

We can see that the sum (solid histogram) of our simulations is in reasonable agreement with 
the ResBos prediction, not only in the prill) dependence but also in the absolute value of the cross 
section. The total cross section is 15.5 pb from ResBos and 13.7 pb from our simulation. The difference 
can be attributed mainly to the lack of finite terms in the NLO correction in our simulation. The 
technical difference in the isolation cut may also have resulted in some difference. 

In our simulation, the contribution from each subprocess is: 4.1 pb (30%) from qq — > 77, 5.1 
pb (37%) from the fragmentation, and 3.7 pb (27%) and 0.8 pb (11%) from the LLL subtracted 
qg — > 77 + q and qq — > 77 + 5 processes, respectively. The contribution from the lowest-order process, 
qq — > 77, is smaller than 1/3 of the total sum. Furthermore, as we can see in Fig. |H1 the prill) 
spectrum of qq — > 77 is apparently different from the sum. The contribution from the LLL subtracted 
qg — > 77 + q is comparable with the lowest-order contribution. However, it must be noted that the 
composition in the above is not physically meaningful. The above is a result when we separate the 
soft and hard radiations at /i = pt of the qq — > 77 or qg — > 7 + q interaction. The composition changes 
if we adopt another definition. In any case, we can at least conclude that the contribution from the 
19 ~ * 77 + 1 process in which the final-state three particles are well separated is significant in the 
non-resonant QED diphoton production. Together with the LL component which is shared between 
the fragmentation process and the LLL subtracted process, non-LL components remaining after the 
subtraction also looks sizable as we can see in Figs. [3] and [H 
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Figure 6: px distribution of the diphoton (77) system in the simulated events satisfying the kinematical 
condition in Eq. ([I]) and the isolation condition. The sum (solid histogram) is presented together with 
the distributions of the subprocesses: qq — > 77, the LLL subtracted qg — > 77 + q and qq — > 77 + q 
and the fragmentation. The ResBos prediction for the same condition is plotted with open circles for 
comparison. 



6 Conclusion 

In conclusion, we have developed an exclusive event generator for non-resonant QED diphoton (77) 
production at hadron collisions, consistently including additional jet productions. The qg — > 77 + q 
process to be included as a radiative process has a final-state QED divergence as well as an initial-state 
QCD divergence. We have developed a subtraction method to regularize the final-state divergence, by 
extending the method developed for initial-state QCD divergences. The subtraction works well as we 
expected. The differential cross section converges to zero after the subtraction at the limit where one 
of the pairs of the final-state 7 and quark becomes collinear and the original cross section diverges. 

We tried to use the PYTHIA simulation for the generation of the fragmentation events, with which 
the subtracted component is to be restored without any double count. We generated 7 + jet events 
by PYTHIA and picked up another photon produced by the parton shower (PS) in the final state. 
We have developed a technique to reconstruct on-shell parton level 77 + jet events from those events 
fully simulated down to the hadron level, in order to test the matching with the events generated 
by GR@PPA. The fragmentation events generated by using the default "old" PS show a reasonable 
matching, though a small mismatch is seen around the boundary and the overall yield seems to be a 
little bit smaller. The treatment of small-Q 2 radiations is still an open question in this simulation. 

A consistent simulation sample was composed by combining the above two generation samples 
with the qq — > 77 and qq — > 77 + g events simultaneously generated with qg — > 77 + q in GRQPPA. 
The event generation was tested for the LHC design condition, pp collisions at 14 TeV, and a typical 
Higgs-boson search condition was required to the events simulated down to the hadron level. A 
typical isolation cut was also applied to the photons. We observed a reasonable suppression of the 
fragmentation events in a collinear region, together with a visible smearing of the boundaries due to 
the application of PS and hadronization. 

The combined events show a behavior which reasonably agrees with the prediction from a re- 
summed NLO calculation by ResBos. The difference in the total cross section at the level of 10% 
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can be attributed to the lack of finite-term contributions in the NLO correction in our simulation. In 
our simulation sample, for which all the energy scales were chosen to be equal to the px of the basic 
qq — > 77 or qg — > 7 + q process, the contribution from the lowest-order qq — > 77 process is smaller 
than 1/3. The contribution from qg — > 77 + q is comparable with it due to a large density of gluons 
inside protons. A consistent inclusion of this process must be necessary for reliable studies. 

A new PS which implements QED radiations together with QCD is under development. We want 
to implement a mechanism to force a hard photon radiation in this PS in order to ensure a reasonable 
generation efficiency for the fragmentation process. 
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